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Abstract: We derive the spectral index of /jvx and its running from isocurvature single 
field and investigate the curvaton models with a negative spectral index of /nl in detail. 
In particular, a numerical study of the axion-type curvaton model is illustrated, and we 
find that the spectral index of Jnl is negative and its absolute value is maximized around 
cr* = 7r//2 for the potential V{a) = m 2 f 2 (l — cos j). The spectral index of Jnl can be 
0{— 0.1) for the axion-type curvaton model. A convincing detection of a positive nf NL 
will rule out the axion-type curvaton model. In addition, we also give a general discussion 
about the detectable parameter space for the curvaton model with a polynomial potential. 
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1. Introduction 



Curvaton model jl], 0, [3| provides an alternative mechanism for generating primordial 
curvature perturbation. One distinguishing feature of curvaton model is that it can produce 
a large local form bispectrum which is proved to be small in the single field inflation model. 
The size of local form bispectrum is measured by the non-Gaussianity parameter Jnl- 
WMAP 7yr data Q implies —10 < Jnl < 74. Recently the curvaton model has been 
widely studied in the literatures §, § 0, |, g, |H| 0, 0, 0, 0, 0, |§ 0, 0, ^ S H 

27], nn, H, II- 



22|, |23|, |2J, |25, 

Even though the curvaton is a light scalar field compared to the Hubble parameter 
during inflation, it still slowly rolls down its potential and its dynamics is expected to 
have some prints on the observational data. Usually the size of local form bispectrum is 
considered to be scale independent. However, it is not a generic prediction for the multi- 
field inflation models, for instance the curvaton model. Recently the authors of 



H|i| 



found that the non-linear evolution of curvaton during inflation can induce a detectable 
scale dependence of Jnl- Such a scale dependence of /jvl is measured by its spectral index 
unl which is defined by 

/ h \ Uf NL 

fNL(k) = f NL (kp) - , (1.1) 



where k p is a pivot scale. Planck 1 35] and CMBPol 36] are able to provide a 1-a uncertainty 



on the spectral index of Jnl for local form bispectrum as follows, [37], 
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where f s k y is the sky fraction. There is no well determined value for the fraction of the 
sky considered for CMB bispectrum measurements. For example, in 0], the authors used 
for the analysis of non-Gaussianity, the KQ75y7 mask, corresponding to a f s k y of 70.6%. 
Roughly speaking, as long as \nf NL ■ /nl\ is not less than 5 or 2.5, the scale dependence of 
Inl w iH be detected by Planck or CMBPol. The large-scale structure data may provide a 
more stringent constraint on the scale dependence of /nl- The scale dependence of J'nl 
will become an important observable in the near future. 

In this paper we mainly focus on the curvaton model and discuss how it can generate 
a negative spectral index of /nl- Our paper is organized as follows. In Sec. 2 we consider 
a general scale dependence of /nl from an isocurvature single field and derive the spectral 
index and its running of /nl- We focus on the axion-type curvaton model where a negative 
rif NL is predicted in Sec. 3. In Sec. 4 we estimate the detectable parameter space for the 
curvaton model with a polynomial potential. Some discussions are included in Sec. 5. 



2. The spectral index and running of f^L from an isocurvature scalar field 

In this section we assume that the curvature perturbation is produced by the quantum 
fluctuation of isocurvature single field a which slowly rolls down its potential at the infla- 



tionary epoch. Based on the so-called 5N formalism [38], the curvature perturbation can 
be expanded to the non-linear order as follows 

C(i/,x) = N a (t f ,ti)6a(ti,x) + -N^{tf,ti)8a 2 {t h yi) + ... , (2.1) 

where N t(T and N )(T(T are the first and second order derivatives of the number of e-folds 
with respect to a respectively. Here tf denotes a final uniform energy density hypersurface 
and ti labels any spatially flat hypersurface after the horizon exit of a given mode. For 
simplicity, U is set to be t*(k) which is determined by k = o(t*)if* for a given mode with 
comoving wavenumber k. Therefore the amplitude of the curvature perturbation generated 
by a takes the form 

P C = K( t *){^)^ ( 2 - 2 ) 
and the non-Gaussianity parameter /tvl is given by 

Here the horizon crossing approximation is adopted. Considering the power spectrum of 
tensor perturbation Pt = 2H 2 /tt 2 , 1 the tensor-to-scalar ratio is defined by 

In order to work out the scale dependence of Jnl, let's introduce a new time t r (> £*) 
which is chosen as a time soon after all the modes of interest exit the horizon during 



x We work on the unit of M p — 1. 
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inflation. If there is a non-flat potential V(a) for a, it will slowly roll down its potential 

3Ha ~ -V'(a), (2.5) 

even though its mass is assumed to be much smaller than the Hubble parameter H. The 
value of a at t r is related to that at time i* by 



t" r da /•*' 

/r, n^y = ~ It, 



dt 



3H{t) ' 



Therefore we have 



(2.6) 

(2.7) 
(2.8) 

(2.9) 

din A; v ' H*dt* 

Here we consider ^F(a r ) = ^^(^§7 + Iff) and 3idcr* = -V'(a*). Taking into 
account that a r is a function of <?*, we have 



<9<7 r 
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One can easily prove 



-F{o T ) 



d 



-F(a r ) = 0. 



da r dN(a r ) 
da* da r 



8 2 a r 8N(ar) ^fda r \ 2 8 2 N(a r ) 



da, 



da 2 



da 2 da r 

Since dN(a r )/da r and d 2 N(a r )/da 2 are scale independent, we find 

d In -/V ,j (t* ) _ 

'/c7(X> 

iV CT (i*) 



din At 
din N :Cra (t* 



— ?7cT<T) 

= 2r/ CT(T + r? 3 - 



dlnA: ' ■" J N <T<T (UY 

where the slow-roll equation for a is adopted and 

= V"(<T*) 

rioo - , 

= V"'(<T*) 

713 ~ 3H 2 ' 

From the above results, the spectral index of and /nl are respectively given by 

dlnP c 



n s = 1 + 



din A; 



= 1 + 2r)aa ~ 2e#, 



and 
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where 

« - -|- (2-18) 

Our result is the same as that in [33 1 . The spectral index of Jnl is proportional to the third 
derivative of potential with respect to a. A free isocurvature single field cannot generate 
the scale dependent Jnl, and the spectral index of /nl is a good quantity to measure the 
self-interaction of such an isocurvature scalar field. 

In addition, the spectral index nf NL may also depend on the scales and its scale- 



dependence is measured by its running which is introduced in [34]. Here we can easily get 



the running of the spectral index nf NL generated by an isocurvature field as follows 

din f 

afNL ^ dhS = {2eH " 11(7(7 ~ m)nfNL " n *^' (2 - 19) 

where 

V'V"" 

* S 31^' (2 - 20) 

Applying the above formula to curvaton model, we get the same result as that in p4| . 
Replacing 2e# by 1 — n s + 2r] aa , the running of spectral index of /nl becomes 

a fNL = (! - n s + Vaa -m~ n f NL ) n fNL- ( 2 - 21 ) 

If rif NL ~ 0(1), otf NL is roughly the same order of rif NL , and then rif NL is not a reliable 
quantity to measure the scale dependence of Jnl any more. In this paper, we mainly focus 
on rif NL and don't plan to extend the discussion about aj NL in detail. 

3. Axion-type curvaton model 

The Nambu-Goldstone boson - axion - with shift symmetry a — > a + b is supposed to be 
a good candidate of curvaton (2|, [|, ||, 19, 2C[. Its potential goes like 



V{a) = mTU " cos - j , (3.1) 
where / is axion decay constant. If a <C /, the potential is simplified to be 

V{a) ~ \™ 2 o 2 + ■•■ ■ (3-2) 

It is roughly quadratic and the spectral index of /nl is expected to be small. However, 
the potential is highly non-quadratic for a ~ / and a large nf NL is expected. 



From the potential in (3T), we get 



where 

9 

m 

(3.4) 
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which is related to 770-0- by 

Vau = r\ mm cos(<r*//). (3.5) 

We see that j] aa = when cr* = nf /2. In order to calculate rif NL , we need to know N )(r (t*) 
and N t<T(T (t*) as well. It is very difficult to get the analytic formula for these two quantities 
in the axion-type curvaton model and we will use the numerical method in 2C ] to compute 



them. The numerical calculation is based on the so-called 5N formalism |J8| where we start 
from any initial flat slice at time U n i, and then the curvature perturbation will be given 
by C(^ x ) = -^(£; x ) — No(t) on the uniform density slicing. Here No(t) is the unperturbed 
amount of expansion. 

After inflation, the vacuum energy which governs the dynamics of inflation was con- 
verted into radiation and the equations of motion of the universe become 

H 2 = \{p r + P„), (3.6) 
Pr + AHp T = 0, (3.7) 

1 

— ( 

2 



Pa = \o 2 + V(a), (3.8) 



a + 3Ha + ^ = 0, (3.9) 
da 

where p r and p a are the energy densities of radiation and curvaton respectively. It is 
convenient to introduce some new dimensionless quantities, such as N(t) = lna(t), x = rat 
and 

6 = a/f. (3.10) 



Now the Hubble parameter becomes 



dN 

H = m—, (3.11) 

dx 



and the above differential equations of motion are simplified to be 



(3.12) 



dN 
dx 

d 2 e | ^dNde | dvje) _ Q 

dx 2 dx dx dO 

where 

V(a) = m 2 f 2 V(6), (3.14) 

and 

pr,ini / Q -. c \ 

a = (3.15) 

and p r ,ini is the radiation energy density at t = ti n i. For the axion-type curvaton model, 

V(9) = l-cos9. (3.16) 

The scale factor can be rescaled to satisfy a(tj n j) = 1 and then N{ti n i) = 0. If the inflaton 
energy decays into radiation very fast, a ~ H 2 n j/m 2 which should be much larger than one. 
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Our numerical computation indicates that the numerical result is insensitive to the choice 
of a for a > 10. In this paper we set a = 10. We also adopt sudden decay approximation 
which says that the curvaton is supposed to suddenly decay into radiation at the time of 
H = r, or equivalently dN/dx = T/m, where T is the curvaton decay rate. This condition 
is used to stop the evolution of curvaton in the numerical calculation. 

In this section we only plan to illustrate the physics for the axion-type curvaton model 
and adopt / = 10 _1 and T/m = 10~ 2 in order to shorten the time for numerical calculation. 
The numerical results of /jvl an d Tif NL are showed in Fig. |J, A similar plot to the bottom 
one in Fig. [l] is obtained for other choice of the parameters in the model. We see that the 
axion-type curvaton model predicts a negative n/ NL . Requiring that the curvaton slowly 
rolls down its potential during inflation implies \rj aa \ <C 1, but r] mm can be ©(HT 1 ~ 1) 
around 6>* ~ ir/2, and the absolute value of nf NL is maximized around 9* = n/2. So nj NL 
can be 0{— 10 _1 ) for axion-type curvaton model. 

If = 7r//2, there is a consistent relation between nj NL and Q-f NL , from Eqs. ( 2.1 S| ) 
and (p|), 



a fNL = (1 - n a - n fNL )n fNL . (3.17) 
The WMAP 7yr data implies n s = 0.963 @, and then a fNL = (0.037 - n fNL )n fNL . 

4. The detectable parameter space for the curvaton model with a polyno- 
mial potential 

In this section we consider the curvaton model with a polynomial potential 

V (a) = lm 2 a 2 + \m A (-Y , (4.1) 



2 \m/ 

where the dimensionless coupling A can be positive or negative. For the axion-type curvaton 
model with a / in Sec. 3, A = — m 2 /(24/ 2 ). The curvature perturbation for this kind 
of curvaton model is discussed in |l3|, [l5|, [jl], ^2], 24, 26, 34] in detail. In particular, the 



authors of [34| only focused on the curvaton model with near quadratic potential and a 



positive dimensionless coupling A. However, in this paper the self-interaction term can be 



subdominant or dominant. Similar to [13], we introduce a new parameter 



(n \ n—z 
— ) (4-2) 

to measure the size of the curvaton self-interaction compared to its mass term. For the 
model with negative s, V (<r*) is required to be positive, or equivalently s > — 2/n; otherwise 
the curvaton field will run away, not oscillate around a = 0. The parameter r} aa is related 
to s by 

n(n — 1) 



m 

V™ = 



3H? 



1 + s 



(4.3) 



Usually en is much smaller than one, and hence a red tilted power spectrum of curvature 
perturbation cannot be naturally achieved in the curvaton model with quadratic potential 



-6- 



2000 
1000 
500 

200 
100 
50 

20 



0.1 



0.2 



0.3 0.4 

ejn 



0.5 



0.6 



0.7 



0.8 



0.6 



* 0.4 



a 
i 



0.2 



0.0 




Figure 1: Jnl and rif NL for axion-type curvaton model. Here we adopt / = 10 1 and T/m = 10 



p^| . Once we take the self-interaction of curvaton into account, this problem is released 
significantly. For s £ [— — , ~ n ( w 2 -i) ]i r] aa < and it can make red-tilted. 
The spectral index of /jvx i n this model can be rewritten by 



5 m 



nfNL GN a f NL > 



(4.4) 
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where 

n(n — l)(n — 2) m? 1 
6 Hi cr* 

for the curvaton model with polynomial potential. Combining with the normalization of 
the curvature perturbation fl2.2|), we obtain 



. 5n(n- l)(n-2) m 2 H* s 



where = yPc, — 4.96 x 10 [[§. In [34], the authors did not consider the constraint 



from the WMAP normalization. We need to stress that the above formula is applicable 
for any value of s > —2/n. Here are three model-dependent free dimensionless parameters: 
m/H*, H*/a and s. Focusing on the curvaton model with near quadratic potential, namely 
\s\ <C 1, N >a is positive and the spectral index of /jvl is negative when s < for a positive 
Jnl, such as that in the axion-type curvaton model. For a curvaton model with n = 4, 
Vmm ~ 10 -2 , H*/a* ~ 1CP 1 and \s\ ~ 1CT 1 , \nf NL ■ /jvx| ~ 6.4 which is detectable for 
Planck. 

Here we also want to estimate the typical value of curvaton field during inflation. In 
a long-lived inflationary universe, the long wavelength modes of the quantum fluctuation 



of a light scalar field may play a crucial role in its behavior [39|, because its Compton 
wavelength is large compared to the Hubble size during inflation. The quantum fluctuation 
can be taken clS cl random walk: 

(« 2 ) = £|t- (4.7) 

On the other hand, the long wavelength modes of the light scalar field are in the slow-roll 
regime and its dynamics is governed by 

3H^ = -m 2 a -nXm 3 (-Y 1 . (4.8) 
at \m/ 

Combining these two considerations, we have 

d(a 2 ) _H^_ 2m 2 2 _ 2nAm 4 (a 2 )"/ 2 



See more discussions in [16, |39|. Typically the vacuum expectation value of a can be 
estimated as 



(7 



VV). (4.10) 



It is reasonable to suppose that cr* sits at the point where the solution of the differential 
equation ( |4.9|) approaches a constant equilibrium value: 



a * = V 87 r 2 (i% S )^' (4U) 



-8- 



with s > —2/n. Now Eq. ( |4.6| ) becomes 

n fNL ■ f N L = sign(N a ) ■ 2.3 x 10 3 n(n - l)(n - 2)(-^) 3 s^l + ^s. (4.12) 

For the typical choice of the vacuum expectation value of curvaton during inflation, the 
signal about the scale dependence of fNL depends on the power of a in the self- interaction 
term, m/H* and s. The detectable parameter space by Planck (\nf NL ■ Jnl\ ^ 5) for 
different n is illustrated in Fig. |2[ The scale dependence of /jvx is detectable if the curvaton 
mass is not too small compared to the Hubble parameter during inflation. If m/ H* < 0.04, 
the curvaton self-interaction term must be positive and dominant in order to obtain a 
detectable scale dependence of fNL- 

5. Discussions 

In this paper we discuss the scale dependence of Jnl from isocurvature single field and the 
curvature perturbation can be expanded as 

Ck = Ct + \fNL{k P )(l + n fNL In ^)(C L * ( L ) k + ... , (5.1) 

where k p is a pivot scale and * denotes a convolution 

(C L * C L )k = J ^C L (q)C L (k ~ q). (5.2) 

The non-linear evolution of isocurvature field induces the scale dependence of /nl- Similar 
to the power spectrum of scalar perturbation, one also introduces the so-called spectral 
index of Jnl to measure the tilt of fNL- We derive not only the spectral index of fNL-, but 
also it running af NL generated by an isocurvature field. The spectral index rif NL is a good 
quantity only when it is small compared to one. 

Focusing on the curvaton model, we find that the spectral index of fNL can be negative 
in the axion-type curvaton model (see Fig. [ij) and in the curvaton model with small negative 
self-interaction term. Even though Fig. |l] depends on the choice of the parameters, such as 
/ and T/m, in the model, we try some different choices and we find that the similar figures 
are obtained. If a positive nf NL is supported by the forthcoming cosmological observations, 
the axion-type curvaton model will be rule out. The spectral index of fNL is a powerful 
discriminate to distinguishing different curvaton models. 

In addition, we also figure out the parameter space for the polynomial potential cur- 
vaton model in which such a scale dependence can be detected in the near future. In our 
consideration, we introduce three parameters, the curvaton mass m, the power of curvaton 
field in the self-interaction term and the relative strength s of self-interaction term com- 
pared to the mass term, to characterize the form of curvaton potential. Combining with 
the estimation of typical vacuum expectation value of curvaton field at inflationary epoch 
and WMAP normalization, we reach Fig. || where a large parameter space for detectable 
scale dependence of fNL is illustrated. We need to stress that our result is applicable even 
when the self- interaction term becomes dominant. 
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Figure 2: The shaded regions illustrate the detectable parameter space for the curvaton model 
with a polynomial potential. 



Finally we also want to point out that the scale-dependent /tvl is not a distinguishing 
feature for the curvaton model. Similar phenomenology can happen in the model [40] 
where the entropy fluctuations are converted into the adiabatic perturbation at the end of 
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multi-field inflation due to the geometry of the hypersurface for inflation to end if there 
are self-interactions for the degrees of freedom along the entropy directions. 
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